Momentum space evolution of chiral three-nucleon forces 
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A framework to evolve three-nucleon (3N) forces in a plane-wave basis with the Similarity Renor- 
malization Group (SRG) is presented and applied to consistent interactions derived from chiral 
effective field theory at next- to-next-to- leading order (N^LO). We demonstrate the unitarity of the 
SRG transformation, show the decoupling of low and high momenta, and present the first inves- 
tigation of universality in chiral 3N forces at low resolution scales. The momentum-space-evolved 
3N forces are consistent and can be directly combined with the standard SRG-evolved two-nucleon 
(NN) interactions for ab-initio calculations of nuclear structure and reactions. 
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The interplay of chiral effective field theory (EFT) and 
Renormalization Group (RG) methods offers new oppor- 
tunities for efficient and simplified microscopic calcula- 
tions of many-body systems based on systematically de- 
rived nuclear forces [IJI2]- Recent applications range from 
calculations of finite nuclei [3H7] to infinite nucleonic sys- 
tems [SI 12] to astrophysical applications 10]. A key chal- 
lenge in all of these cases is the control of many-body 
forces. In this letter we present a new tool to address this 
challenge, the RG evolution of 3N forces in a plane-wave 
basis (see Fig. [T]), which we use to explore universality in 
the evolved chiral 3N forces (see Fig. |4]). 

The SRG provides a framework to construct unitary 
transformations that consistently renormalize all opera- 
tors, including many-body forces, while preserving low- 
energy observables [II1II2]- Wegner's formulation of the 
SRG is a continuous series of unitary transformations of 
the Hamiltonian H = Tioi + as a function of the flow 
parameter s: 



ds 



(1) 



Here rjg specifies the unitary transformation, Tici denotes 
the relative kinetic energy and V all the interparticle in- 
teractions. In the following we make the common choice 
Vs — [Trci, Hs], which generates transformations that lead 
to a decoupling of low and high momenta in NN interac- 
tions as they are evolved to lower resolution scales [13]. 
Such an evolution leads to much less correlated wave 
functions at low resolution and the nuclear many-body 
problem becomes more perturbative. 

In recent RG-based calculations, two different strate- 
gies have been used to handle the 3N forces. Starting 
from nuclear NN and 3N forces, derived and fitted in 
chiral EFT, it is possible to systematically evolve the full 
Hamiltonian. So far, this has been achieved by represent- 
ing Eq. ^ using a discrete harmonic oscillator basis [13] ■ 
Results for light nuclei based on such evolved interactions 
are promising For heavier nuclei however, significant 
scale dependencies have been found [U [15] , which suggest 
that infinite matter will not be realistic. These could be 
indications of significant induced 4N forces or of an in- 
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FIG. 1. (Color online) Ground state energy of as a func- 
tion of the flow parameter A = s~^^'^ for different chiral in- 
teractions at N'^LO, which are labeled by the value of the 
cutoffs A/ A (see Ref. [19]). The solid lines show the results 
for A/A = 550/600 MeV separately for NN-only, NN-(-3N- 
induced and NN-f3N-full (see Ref. [H]). The dashed fines 
show the NN-|-3N-full results for the other given values of 
A/A. Only np NN interactions have been used. 



sufficient evolution of 3N forces due to basis truncations. 

An alternative strategy has been used for calculations 
of infinite nuclear systems [H] (and elsewhere [51|S]). In- 
stead of fitting the 3N forces at the chiral EFT cutoff 
scale, only the NN interactions have been evolved with 
RG methods and then the short-range parameters of the 
N^LO 3N forces have been determined from fits to few- 
body systems at the low-momentum scale. This proce- 
dure assumes that the long-range part of the 3N forces 
remains invariant under the RG transformations and that 
the N^LO operator structure is a sufficiently complete op- 
erator basis that induced contributions can be absorbed 



2 



to good approximation. The results from this procedure 
are found to be in agreement with nuclear phenomenol- 
ogy within the theoretical uncertainties [8l[9j. However, 
these results do not imply that induced 3N forces have 
the same structure as chiral N^LO 3N interactions. 

A complementary evolution of nuclear many-body 
forces using the SRG in a momentum plane-wave basis 
would allow us to reconcile the seemingly inconsistent re- 
sults of these strategies. In this letter, we report the first 
results for the evolution with Eq. ([T]) of realistic nuclear 
NN and 3N forces in such a basis. A similar approach 
has been used recently for the evolution of bosonic inter- 
action models in one dimension |16j . As we will discuss 
below, however, the coupling of partial waves makes the 
RG evolution more challenging in three dimensions. 

We adopt the notation and conventions of Ref. [T3| and 
write the interaction in the form Vg = Vu + V13 -I- V23 + 
V123 and the kinetic energy as T^d = T12+T3 = T13+T2 = 
T23+T1. Here Vij denotes NN interactions between parti- 
cle i and j and V123 the irreducible three-body potential. 
The kinetic energy terms Tjk and Ti correspond to the 
contributions of the Jacobi momenta pi and 9^, respec- 
tively. We recast Eq. ([T]) as separate RG flow equations 
for the two- and three-body interactions [13]: 



ds 

dVl23 

ds 



[[Tij, Vij] , Tij 



(2) 



[[^12, ^12], 1^3 + "1^23 + ^^123] 

+ [[Tu,Vu],Vi2 + V23 + Vi23] 
+ [[^23, ^^23] , V"l2 + Vi3 + V123] 
+ mohVl23],Hs]. 



(3) 



Compared to Eq. Q, this system of differential equa- 
tions has the important advantage that terms result- 
ing from spectator particles in two-body interaction pro- 
cesses have been eliminated explicitly. In a momentum 
basis, the spectator particles lead to delta functions that 
make the representation of Eq. ([T]) problematic [T7] . 

The flow equation ^ is solved in a two-body basis and 
then embedded in Eq. ^ by using 



132V23-P132, 



-P123 ^23^123 ■ 



(4) 



Here P123 = ^12^23 (^132 = P13P23) permutes three par- 
ticles cyclically (anti-cyclically) with Pij denoting two- 
particle transpositions. We represent Eq. ([2]) in a stan- 
dard partial- wave basis of the form \p; {LS)JT), where 
L, S, J and T denote the orbital angular momentum, 
spin, total angular momentum and isospin of the inter- 
acting pair of particles with relative momentum p. For 
the three-body basis we choose [TTl [18] 



(5) 



\pqa), = \p,q,; [{LS)J{ls,)j] JMTt,)T%) 



where pi and qi denote the three-body Jacobi momenta 
of particle i. The quantum numbers I, Si = 1/2, j 



and ti = 1/2 label the orbital angular momentum, spin, 
total angular momentum and isospin of particle i rela- 
tive to the center-of-mass of the pair with momentum p. 
J and T denote the total three-body angular momen- 
tum and isospin, which are equal to 1/2 in the three- 
nucleon bound states. For details we refer the reader to 
Refs. [T71 HH]. The basis states ([S]) are not completely 
antisymmetric. It is most natural to evolve the antisym- 
metrized interaction 

(M«|Fi23b'g'a') = .(M«|-4i23V^/23^123|p'g'a')» ' (6) 

with ^123 = (1 + -P123 + -P132) and y/23 being the i-th 
Faddeev component of the three-body interaction |17j . 

We solve Eqs. ([2]) and ([3]) for a set of five different ini- 
tial NN and 3N interaction potentials derived from chiral 
EFT at next-to- next-to-leading order (N^LO) [IS]. The 
3N forces low energy constants cd and ce have been fit- 
ted for five different cutoff combinations A and A to the 
•^H binding energy and the nd doublet scattering length 
^a„p [ini HO]. For this first exploratory work we rep- 
resent the initial Hamiltonian at s = in a basis with 
— 26 three-body partial- wave channels ( J < 3). The 
corresponding ground state energies are shown in Fig. [T] 
They are found to be within 12 keV of the converged val- 
ues based on the antisymmetrized interaction (|6]). We 
neglect the effects of charge symmetry breaking, i.e., we 
use the neutron-proton interaction in all isospin channels 
of the NN force. Due to this we obtain ground state 




FIG. 2. (Color online) Ground state energy of '^H as a func- 
tion of the flow parameter A = s^^''* for different SRG model 
space sizes and a fixed model space for the initial Hamil- 
tonian A^^ — 5 (for a detailed definition of the 3N channels 
see Ref. [U])- The dotted line shows the energy at A = 00. 
Only np NN interactions have been used. 
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FIG. 3. (Color online) Contour plots of the evolved 3N potential (^a = l|V^i23iC'Q' = l) Bit the hyperangle 9 = -k/12 (see 
Ref. [17] for definition of partial waves). The upper panel shows the potential with A/A — 550/600 MeV and the lower panel 
A/A = 450/500 MeV. 



energies that are about 150-260 keV more bound than 
the experimental value Egs = -8.481821(5) [H] for all 
studied potentials. For future applications it will be pos- 
sible to incorporate the charge dependence of the nuclear 
interactions, as well the Coulomb potential, in the flow 
equations. 

The SRG evolution is performed in a larger model 
space with Na = 42 > channels (J < 5). This basis 
size is also used for the solutions of the Faddeev equations 
at the different resolution scales. As can be seen in Fig.jl] 
the ground state energy is not perfectly invariant under 
the RG transformations. This is because a finite basis of 
the form ^ is not complete under cyclic and anticyclic 
permutations of particles, as the permuation operator 
i{^pqa\Pijk\p'q'a'')i couples in general all partial waves 
a and a' [T71[T5]. As a consequence, non- vanishing ma- 
trix elements in all three-body partial waves are induced 
when two-body operators are embedded in a three-body 
momentum basis ([5]) via Eq. (|4]) . This problem is absent 
in one dimension ^16j or in a discrete oscillator basis |14j , 
where the permutation operator is block-diagonal in a 
given model space of size A'max- The violation of uni- 
tarity is correlated with the value of the cutoff values 
A/A. For intermediate values A/A = 550/600 MeV we 
find a maximal variation of 5 keV of the ground state en- 
ergy, compared to about 400 keV if induced 3N forces 
are completely neglected. For small values A — 450 
MeV the variation is even smaller (< 2keV), whereas 
for A = 600 MeV it is about 16 keV. This is natural since 
high momentum modes couple stronger with higher par- 
tial waves under permutation of particles. A signifcant 
part of the variations for large A results from an insuffi- 



cient antisymmetrization of the initial 3N forces based on 
a partial- wave representation of A123 (see Eq. For 
future applications it will be possible to optimize the evo- 
lution by using exactly antisymmetrized initial 3N forces 
(see Refs. 

The RG evolution can be improved systematically by 
increasing the model space size Na ■ This is demonstrated 
in Fig. [2j For illustrative purposes we restrict the initial 
Hamiltonian at s = to a small model space with iV° = 5 
channels, which includes the ^Sq and '^Si-'^Di pair partial 
waves (see Ref. QJj). The figure shows the triton ground 
state energy for different SRG model space sizes Na as 
a function of the resolution. For the largest basis size 
shown, the resolution dependence is smaller than 1 keV, 
in contrast to about 30 kcV for the smallest basis. 

Low-resolution NN interactions have been found to be 
quantitatively very similar [21 [13] . This universality can 
be attributed to common long-range pion physics and 
phase-shift equivalence of all potentials, which is reflected 
in the matrix elements at low resolution. It is an interest- 
ing question if the same is true for 3N forces since there 
are important differences: First, 3N forces are fixed by 
fitting only two low-energy constants cd and c^;, in con- 
trast to numerous couplings in NN interactions. Second, 
3N forces give only subleading contributions to observ- 
ables. Since universality is only approximate in NN in- 
teractions, it is not obvious to what extent 3N forces are 
constrained by long-range physics at low resolution. 

In Fig. [3] we show representative examples of 3N forces 
matrix elements at different resolution scales. To do so 
we introduce the hyperradius = + and the hy- 

perangle tan0 = 2p/{\/3q) and choose an arbitrary value 
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FIG. 4. (Color online) Matrix elements of the initial 3N forces (dotted lines) compared to evolved 3N forces (solid lines) at 
A — 1.5 fm~^ for different interactions labeled by the values of the cutoffs A/A. We choose p' = 1.0 fm~^, q' — 1.25 fm~^ and 
a — a' — 2 (see Ref. [I7]). In the left panel we fixed q = 1.5 fm~^ and in the right panel p — 0.75 fm~^. The shaded band 
marks the maximal variation between different evolved 3N force matrix elements. 



6 — 7r/12. The upper panel shows the matrix elements of 
the channel a = a' = 1 for A/A 550/600 MeV and the 
lower panel for A/A = 450/500 MeV. At large resolution 
scales A = cx) (s = 0) the off-diagonal couplings are much 
stronger for larger cutoffs A and the potential is in general 
more repulsive. As we evolve to lower resolution, the re- 
pulsive off-diagonal couplings get successively suppressed 
and finally at A = 1.5 fm~^ we find matrix elements of 
significant size only below < 2fm~^ and around the 
diagonal ^ ~ These features are general and hold 
for matrix elements at all hyperangles and partial waves. 
In addition, the overall effects of the SRG evolution are 
stronger for initial potentials with large cutoffs A/A. In 
summary, the evolved 3N matrix elements show the same 
tendencies found in evolved matrix elements of NN inter- 
actions [13]. 

In Fig. [4] we show in more detail initial and low- 
resolution (A = 1.5 fm""'^) matrix elements for a = a' ~ 2 
(see Ref. fTT) at some typical fixed momenta p\q' ^ 
Ifm^^ for all five different chiral interactions. We find 
a remarkably reduced model dependence for evolved in- 
teractions in this kinematical region. We find that this 
approximate collapse of matrix elements is most pro- 
nounced in the 3N channels in which the three-body con- 
tact interaction ce contributes {a, a' = {1,2}). This 
suggests that the coupling constant c_e; is flowing to an 
approximately universal value at low resolution. In addi- 
tion, new momentum-dependent universal structures are 
induced at low resolution, as can be seen in the right 
panel. If one or more momenta becomes very small the 
model dependence of the matrix elements tends to be- 
come larger. However, the phase space of these matrix 
elements is suppressed and are less relevant for physical 
observables. We emphasize that all these observations 



are based on initial potentials that have been derived at 
N^LO in chiral EFT. At this order, phase shifts at higher 
energies are not as well described as at N^LO, which re- 
duces the degree of universality of NN interactions at 
low resolution. It will be very interesting to investigate 
if universality of 3N forces becomes more pronounced by 
including contributions at N'^LO PHHSj. 

The low-resolution interactions derived here are ideal 
for microscopic calculations of nuclear systems. Evolu- 
tion in a momentum basis has several advantages com- 
pared to a discrete oscillator basis. First, the oscillator 
basis has intrinsic infrared and ultraviolet cutoffs that 
depend on the basis size and oscillator parameter hfl ^ , 
which lead to convergence issues for 3N forces. The 
problems can be avoided by first evolving in momentum 
space and then using a straightforward transformation 
to an oscillator basis with any H^l. This enables new 
tests of SRG interactions in finite nuclei within the no- 
core-shell model [55] and coupled cluster [37] that al- 
low systematic studies of the role of induced 4N forces. 
Second, the momentum-space interactions can be used 
directly in calculations of infinite systems within many- 
body perturbation theory. This will test whether con- 
sistently evolved NN plus 3N forces, initially fit only to 
few-body properties, predict empirical nuclear saturation 
properties within theoretical errors, as found previously 
for evolved NN forces combined with fitted 3N forces [3] . 

Finally, since SRG transformations are usually charac- 
terized by the coupling of momentum eigenstates, the 
momentum basis is a natural basis in which to con- 
struct the SRG generator 77^. In particular, momentum- 
diagonal generators such as T^^i (as chosen here) can be 
implemented very efficiently in a momentum basis and 
it is straightforward to generalize to the Hamiltonian- 



5 



diagonal form advocated by Wegner [T^]. Other gener- 
ators can significantly improve the efficiency of the RG 
evolution [28] or achieve alternative RG decoupling pat- 
terns, such as a flow towards a block-diagonal Hamilto- 
nian [25] . The possibility of using the generator to sup- 
press the growth of many-body forces is also under active 
investigation. 
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